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Abstract 

We study black hole solutions in asymptotically ^^5*4 spacetime with scalar hair. Fol- 
lowing AdS/CFT dictionary these black holes can be interpreted as thermal states of 
2+1 dimensional conformal gauge theory plasma, deformed by a relevant operator. We 
discover a rich phase structure of the solutions. Surprisingly, we find thermodynam- 
ically stable phases with spontaneously broken global symmetries that exist only at 
high temperatures. These phases are metastable, and join the stable symmetric phase 
via a mean-field second-order phase transition. 
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1 Introduction 



The power of holographic gauge theory/string theory correspondence of Maldacena [1] 
is that it extends beyond the conformal examples, and thus can be used as a tool to 
study critical phenomena in strongly coupled systems. Furthermore, while there are 
few explicit realizations of this duality in a non-conformal setting, the application of 
the effective Maldacena framework provides a new valuable (phenomenological) tool in 
addressing the strongly coupled dynamics. Such a framework was applied with intrigu- 
ing results to formulate the holographic model of QCD [2], non-relativistic AdS/CFT 
correspondence [3, 4], effective description of viscous hydrodynamics [5, 6], holographic 
theories of superfluidity [7] and superconductivity [8], to name some. Often a posteri- 
ori, these phenomenological models of the Maldacena correspondence were embedded 
in string theory. The application of the holographic framework itself, rather than its 
specific string theory realization, is particularly relevant to the physics of the continu- 
ous critical phenomena where one expects organization of the physical phenomena in 
universality classes, thus obliterating the importance of any particular representative 
in a universality class. While it is certainly interesting to find embedding of the models 
discussed here in string theory, the latter is not our primary objective in this paper. 

In this paper we follow up on a very interesting suggestion by Gubser^ [9] and 
will attempt to engineer finite temperature critical phenomena in non-conformal 2+1 
dimensional strongly coupled systems^. We now proceed to discuss the motivation of 
our model from the field-theoretic perspective. Consider a 2+1 dimensional relativistic 
conformal field theory. A deformation of this CFT by a relevant operator Or, 

Ti-CFT — H — HcFT + KOr (l.l) 

softly breaks the scale invariance and induces the renormalization group fiow. Our 
experience with similar deformations [12, 13] suggests that in a thermal state Or will 
develop a vacuum expectation value. This expectation value tends to zero in a high 
temperature phase, but it grows as the temperature T becomes comparable to a scale 
A, set by the coupling A^. We assume that the deformed theory Ti has an irrelevant 
operator O,. Note that we do not (further) deform 7i by Oi, as this would destroy a 

^We would like to thank Omid Saremi for bringing this work to our attention and for valuable 

discussions concerning his related work [11]. 

^Our choice of the dimensionality was dictated by the further possibility to study these systems 
(in a holographic setting) in external magnetic fields [10]. 
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UV fixed point defined by TicpT- We consider the case in which the operator Oi has 
an associated Z2 discrete symmetry. As long as the operators Or and do not mix 
under the RG flow, the expectation value of Oi is zero. Of course, the generic situation 
is that the operator mixing is present, in which case Oi will develop a vev in a thermal 
state. The absence of discrete symmetries would typically produce a nonvanishing 
expectation value of irrelevant operators even in the high temperature phase [14, 15]. 
On the other hand, it is natural to expect that a discrete symmetry (preserved by a RG 
flow) would be protected at least when the operator mixing is small, i.e., in the high 
temperature phase. As T < A one expects that the mixing would be large, producing 
the condensate of Oi and thus spontaneously breaking Z2 symmetry. A holographic 
implementation of this phenomenon was proposed in [9] . 

Rather surprisingly, wc find that in our specific holographic realization of Gubser's 
proposal [9] , spontaneous breaking of a discrete symmetry occurs at high temperature. 
The thermal phases with the broken symmetry exist only above certain critical tem- 
peratures — different for each distinct broken phase. In a holographic setting, each 
broken phase is identified by the number of nodes in the wavefunction of the super- 
gravity field dual to an irrelevant operator Oi. Although thermodynamically stable, 
these phases have higher free energy than the symmetric phase. Some of them connect 
to a symmetric phase via a second-order mean-field transition, with {Oi) as an order 
parameter. 

In the next section we discuss our holographic model in detail. The results of the 
thermodynamic analysis are presented in section 3. We conclude in section 4 with 
future directions and further conjectures. 



2 A holographic model of an exotic critical phenomena 

In this section wc discuss an explicit realization of Gubser's holographic model of 
critical phenomena [9]. 

The effective four-dimensional gravitational bulk action, dual to a field-theoretic 
setup discussed in the introduction, takes form^ 

5-4 ^ScFT + ^r + Si = ^J dx^^ [CcFT + J^r + A] , (2-1) 

jCcft = + 6 , A = (V^)' + 0^ A = (Vx)' - 2x' - g(l>\' , (2.2) 

^We set the radius of an asymptotic ^1^5*4 geometry to unity. 
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where we split the action into (a holographic dual to) a CFT part Soft] its deformation 

by a relevant operator Or', and a sector involving an irrelevant operator Oi along 

with its mixing with O,- under the RG dynamics. The four dimensional gravitational 

constant n is related to a central charge c of the UV fixed point as 

192 , , 

c= — . (2.3) 

In our case the scaling dimension of Or is either 1 or 2, depending on which of the 
two normalizablc boundary modes of a gravitational scalar [16] we treat as a field- 
theoretic parameter (and thus keep it constant as we vary the temperature), and the 
scaling dimension of Oi is 4. In order to have asymptotically AdS^ solutions, we 
assume that only the normalizablc mode of Oi is nonzero near the boundary. Finally, 
we assume that g < in order to holographically induce the critical behavior [9]. 

Effective action (2.1) has a Z2 x Z2 discrete symmetry that acts as a parity transfor- 
mation on the scalar fields and x- The discrete symmetry — > — is softly broken by 
a relevant deformation of the AdS4^ CFT; while as we will see, the x ~^ ~X symmetry 
is broken spontaneously. 

We will study Schwarzschild black hole solutions in (2.1) with translationary in- 
variant horizons and nontrivial scalar hair. Thus, we assume: 

dsi^-ci(ry dt'' + C2(r)^ [dxl + dxl]+cldr\ = 0(r) , x = x(r) • (2.4) 

As in [14] , for the precision numerical analysis we find it convenient to introduce a new 
radial coordinate x as follows 

1-x^^, (2.5) 

so that a; — > corresponds to the boundary asymptotic, and y = 1 — a; — > corresponds 
to a regular Schwarzschild horizon asymptotic. Further introducing 

the equations of motion obtained from (2.1), (2.2) with the background ansatz (2.4), 
and the boundary conditions discussed above, imply 



a —a 



--pi x'/' + p2 + —pIx + C»(x^/3) , (2.7) 
20 



near the boundary x — > 0+, and 

a = a (ao" + + 0{y')) , <P = p'o + 0{y') , cj> = ^ + 0{y') , (2.8) 

near the horizon y — 1 — x ^ 0+. Apart from the overall scaling factor a (which is 
related to the temperature), the background is uniquely specified with 3 UV coeffi- 
cients {pi,P2,X4} and 4 IR coefficients {cq, aJ^jPg, Cq}. Recall that both modes of 
near the boundary are normalizable. Thus we have a freedom to interpret the relevant 
deformation (1.1) as due to Or with dim[Cr] = 2, in which case we have to keep the 
combination pia = A fixed, or due to Or with dim[Cr] = 1, in which case we have 
to keep = A^ fixed. We can determine these scalings by carefully matching the 
boundary asymptotic finite temperature background (2.4) to the boundary asymptotic 
of the zero temperature background — in the latter case it is clear what is meant to 
keep Xr fixed^. As we point out below, the highly nontrivial check on our identification 
is the consistency of the hairy black hole thermodynamics with the first law of ther- 
modynamics. The scale A is the scale set by the relevant couphng Xr in (1.1). Notice 
that the two solutions 

(/>i,P2,X4,ao,a^>Po,Co) & (-Pi, -P2, -X4, Oq, a?, -Po, -Co) 

are equivalent. Thus we can assume > for dim[(9r] = 2, andp2 > for dim[(9r] = 1- 
It is straightforward to compute the temperature T and the entropy density s of 
the black hole solution (2.4) 

SttTV _ 6(a{i)^(6 - 2(A;i)^ + (pi^f Mfj^m ^ ^3.9) 



a ) 2>a1 -I- Oq 

s = ^^^i:a^ {alf . (2.10) 

The energy density E and the free energy density T can be computed using the tech- 
niques of the holographic renormalization [17] . There is a slight subtlety in holographic 
renormalization in our case associated with the fact that both modes of Or are nor- 
malizable near the boundary [18]: 



^Identical procedure has been used in [14]. 
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when dim [Or] = 2 we have 



(2.11) 



when dim [Or] = 1 we have 



c 



3"^ 2 y 3a^ + 



384 ^ / 1 



(2.12) 



In both cases we independently evaluated {J^,£,s}. A nontrivial check on the holo- 
graphic renormalization is the automatic fulfillment of the basic thermodynamic rela- 
tion 

JT = £: - Ts . 

The UV fixed point is described by AdS4 Schwarzschild-black hole solution 

"^^^ ^ (2x-x^)V3 ' ^(^) = ^(^) = ' (2.13) 
in which case we have 

4^' (^~~27' (^~^' ( ^ 

We conclude this section with comments on the numerical procedure. We use 
numerical technique identical to the one developed in [14]. For dim[(9,.] = 2 we vary 
Pi, for each given value of pi we compute the remaining UV and the IR coefficients 
{P2, X4i ^oi ^iiPo' '^o)- Notice that there are just the right number of these coefficients 
to uniquely determine the solution of a system of 3 second order differential equations 
(the equations of motion) for the background functions {a{x) , (/)(x) , x{^)}- Given the 
numerical data, we can use (2.9) and the condition pia = A to present the data, say 
for the free energy, as ^ versus ^. The procedure is then repeated for dim[(9r] = 2 
with appropriate replacements of pi with p2. 
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Figure 1: The free energy density of hairy black holes in asymptotic AdS4 geometry for 
different deformations of the UV fixed point. The red curves correspond to symmetric 
phases with (Oj) = 0. The other curves correspond to phases with spontaneously 
broken Z2 symmetry, (Oi) 7^ 0. The purple curves represent symmetry broken phases 
with the holographic condensate wavefunction x{^) ^^ot having any nodes for x G [0, 1]. 
The green curves represent symmetry broken phases with the holographic condensate 
wavefunction x(^) having exactly one node for x G [0, 1]. 

3 Thermodynamics of hairy black holes in asymptotic AdS^ 
geometry 

In this section we present results of the thermodynamic analysis of the hairy black 
holes. We choose g = —100 in (2.2). 

Figures. 1 and 2 represent the free energy density and the energy density (cor- 
respondingly) of the hairy black holes in asymptotic AdS^ geometry. Each plot is 
labeled by the dimension of the relevant operator Or used to generate the holographic 
RG flow, see (1.1). The red curves correspond to Z2 symmetric phases with a vanish- 
ing expectation value of the irrelevant operator Oi (recall in our model dim[(9j] = 4 
). The remaining curves represent phases with the broken Z2 symmetry, and the non- 
vanishing expectation value for O^. The purple curves denote symmetry broken phases 
in which the wavefunction of the gravitational holographic dual to Oi — x{^) — has 
no nodes for x G [0, 1]. The green curves denote symmetry broken phases in which the 
wavefunction x{x) has exactly one node for x G [0, 1]. 

Some of the comments below follow from the data presented; some other comments 
rely on the results of the analysis we omit here. 
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Figure 2: The energy density of hairy black holes in asymptotic AdS^ geometry for 
different deformations of the UV fixed point. The color coding is as in figure 1. 

. dim [a,] = 2 2 dim [a] = 1 




Figure 3: The speed of sound of hairy black holes in asymptotic AdS^ geometry for 
different deformations of the UV fixed point. The color coding is as in figure 1. 



■ We verified the consistency of our analysis and the identification of the mass scale 
in deformation (1.1) by checking the first law of thermodynamics: 



dj^ 



-sdT . 



In our numerical results 



1 + 



dJ^ 



sdT 



< 10^ 



(3.i: 



■ All the symmetry broken phases for the deformation with dim[(9r] = 2 (and some 
of the phases generated by the dim[(9r] = 2 RG flow) exist only above certain critical 
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A - In ( ^ -u M 

^ - I {7rT)3 + 27 



dim[0^] = 2 



10 II 12 



ln(-^) 



Figure 4: The difference between the free energy density of the "purple" symmetry 
broken phase and the symmetric phase (see figure 1) for ^ ^ as a function of 
ln(— (yf). The straight hne fit to the data produces the following functional dependence: 
A ^ 2.32861 - 1.01493 ln(-^). 

temperature, Tc. Notice that this is different for phases in which the condensate 
x(x) has different number of nodes — j'vurpie ^ rpgreen ^j-ecaAl that in purple phases 
x{x) does not have any nodes, and in green phases x(^) has one node ). 

■ All the symmetry broken phases have a higher free energy density compare to a 
symmetric phase. Thus, at best, they are metastable. Further work is needed to deter- 
mine if phases with spontaneous symmetry breaking are stable [19]. Note that if these 
phases are shown to be perturbatively unstable (having a tachyon in the spectrum of 
fluctuations), it will be rather unusual given that these phases survive to arbitrary high 
temperatures ^ ^ 0. 

■ All the broken phases for the dim[(9r] = 2 deformation, and the phases for the 
dim[(9,.] = 1 deformation which "end" on the symmetric phase (the red curves), un- 
dergo a mean-field second-order phase transition. In all 

(a)ocx4 0c(T-^e)^/^ (3.2) 

+ {T — TY as T > T 

jrbroken j^symmetric ^ "^z ' ' (3 3) 

as T < Te . 

■ All the phases we discuss are thermo dynamically stable, at least in the temperature 
range under consideration. Figure 3 represents the dependence of the speed of sound 
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squared cj. on ^. 

■ We found symmetry broken phases with x{^) having more than one mode for x G 
[0, 1], thus we beheve there are infinitely many high temperature phases with {Oi) ^ 0. 
We expect that all these phases share the same general features as the phases we 
discussed in details above. 

■ The data in figures 1-3 correspond to RG flows in which the mixing between Or and 
Oi (the coupling constant g in (2.2)) is fixed, and negative. We studied the dependence 
of the thermodynamics on g. As \g\ decreases, the broken phase free energy density 
evaluated at a given temperature increases. This increase is particularly dramatic 
as g ^ —4 ■ ■ ■ — 3 (for a purple phase) suggesting that there is a critical mixing gc 
(presumably different for different broken phases) at which this symmetry breaking 
thermodynamic phase disappears. We did not find any symmetry broken phases for 
positive g, thus it is likely that gc < 0. Perhaps even more interesting is the behavior of 
the symmetry broken phases for large mixing between Or and O^, i.e., when g — oo. 
In figure 4 we show that as the mixing between the operators along the RG flow 
increases, the broken phases approach the symmetric phase at high temperatures. The 
latter fact is evident both in the behavior of the free energy density and in the vanishing 
of {Oi) for the high temperature broken phases as — > — oo. 

4 Conclusion 

In this paper we used the ideas of Gubscr [9] to holographically engineer RG flows in 
strongly coupled 2 + 1 dimensional gauge theory models with exotic flnite temperature 
phases. Speciflcally, we considered a relevant deformation of a relativistic CFT^ with 
an operator Or, dim[(9j.] = {1,2}. We allowed a mixing of an irrelevant operator Oi, 
dim[(9j] = 4, with Or along the RG flow. We discovered a plethora of exotic phases 
in the system at flnite temperature. We found the phases that exist only at high 
temperature, spontaneously break Z2 symmetry of the model and join the symmetric 
phase with a mean-fleld second-order phase transition. All the phases we discuss are 
thermodynamically stable. The phases with the broken symmetry have a higher free 
energy density than the symmetric phase; thus, at best, they are metastable. 

Still a lot of work (both technical and conceptual) is needed to understand these 
exotic phases. 

■ It would be interesting to study perturbative stability of the phases with the bro- 
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ken symmetry. Using gauge theory/ string theory correspondence it is easy to argue 
that thermodynamic instabihties of the gravitational backgrounds reveal themselves as 
(perturbative) dynamical instabilities [20] as well. The reverse statement is not neces- 
sarily true. For such analysis one needs to compute the spectrum of the quasinormal 
modes associated with the fluctuations of 

■ We made a lot of conjectures about the phase structure of the theory based on ex- 
tensive numerical work. It would be interesting to analytically prove the disappearance 
of the symmetry breaking phases for small mixing, g > Qc- Likewise, it is important 
to show analytically that there is no phase crossing in the system as g ^ — oo. If such 
a phase crossing does occur, the high temperature symmetry broken phases would be 
thermodynamically preferable; in addition, there could be interesting critical behavior 
associated with such crossing. 

■ We did not study the question of embedding the RG flows discussed here in string 
theory — it might very well be that there is an (unknown to us at this stage) obstruc- 
tion for such an embedding. It would be interesting to see if similar phenomena occur 
in string theory realizations of holographic dualities of seemingly related systems [21]. 

■ We want to emphasize once again the unusual features of the phases with the broken 
symmetry. It is well known that a two-dimensional Ising model allows for a dual re- 
formulation in terms of a "disorder" operator^. The Hamiltonian of the Ising model is 
local when expressed in terms of the disorder operators. While the order parameter in 
the spin formulation vanishes above certain temperature, in the disorder formulation, 
the disorder parameter (probing the condensation of the kinks) vanishes in the low 
temperature phase, akin to what we see in our model here. The relation between the 
original spin and the dual disorder operator in the Ising model is nonlocal. In our 
model both the operators Or and Oi are local. If Oi is interpreted as an analog of 
the disorder operator in the Ising model, how can it be mutually local with Or — an 
analog of a local operator in the spin formulation of the Ising model? 

■ There are known condensed matter systems which have a disordered high temper- 
ature phase [23]. It would be interesting to study if a close analogy between those 
systems and the holographic models discussed here could be established. 

■ Finally, it would be interesting to develop in more details the idea of 'microengineer- 
ing' the renormalization group flow by specifying not only the operators that deform 
the (UV) fixed point of the theory, or simply are present in the theory, but also impos- 

^See [22] for a review. 
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ing — by hand — the operator mixing under the renormahzation group flow. While 
this is quite unusual in a standard field-theoretic setting, it is a real possibility in the 
holographic constructions. 
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